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Abstract. The Hipparcos mission has provided very high 
quality parallaxes of a sample of a hundred nearby disk 
stars, of spectral types F to K. In parallel, bolometric 
fluxes, effective temperatures, and accurate Fe/H ratios 
of many of these stars became available through infrared 
photometry and detailed spectroscopic analyses. These 
new accurate data allow to build the Hertzsprung-Russell 
diagram of stars of the solar neighbourhood with the 
smallest error bars ever obtained. 

We analyse these observations by means of theoreti- 
cal stellar models, computed with the most recent input 
physics. 

We first examine the positions of the objects versus 
standard theoretical isochrones, corresponding to their 
chemical composition and age. For these isochrones we 
have first assumed that the helium content was varying 
in locksteps with metallicity. The comparison becomes 
age-independent in the lower part of the HR diagram, 
where evolutionary effects are negligible. We show that 
for the unevolved stars, the agreement between real stars 
and models is fairly satisfactory for stars with metallicity 
within ± 0.3 dex of the solar metallicity, but that a conflict 
exists for stars with metallicity less than [Fe/H] = —0.5. 
This conflict cannot be resolved by decreasing the helium 
abundance: values of this abundance below the primordial 
abundance would be required. 

On the basis of recent works, we show that the addition 
of two processes not included in standard models can help 
solving the above discrepancy. These are (i) correcting the 
LTE iron abundances using a non-LTE approach and (ii) 
including microscopic diffusion of He and heavier elements 
in the stellar interior. The case of the binary star /i Gas is 
particularly useful to support this conclusion as its mass 
is also known from its orbit. After inclusion of the two 



effects, fi Gas A falls on its expected isochrone, within the 
error bars corresponding to its mass. 

AU stars with -0.3 < [Fe/H] < 0.3 are located be- 
tween the helium-scaled isochrones corresponding to these 
metallicities. However five of them are not located exactly 
where they are expected to be for their metallicity. This 
may reflect a helium content lower than the metallicity- 
scaled value. But not necessarily, as a possible sedimenta- 
tion of the elements might complicate the determination 
of the helium content . The age of main sequence solar 
composition stars covers a large range, and the effects of 
sedimentation are time dependent. 
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1. Introduction 

Valuable tests of the stellar evolution theory and strong 
constraints on the physical description of stellar interiors 
are mainly provided by the study of the most-accurately 
observed stellar objects. The number of observable pa- 
rameters accessible through observations is important as 
well as the accuracy on their determination. The Sun, the 
nearby stars including members of visual binary systems, 
and stars belonging to open clusters represent some of 
these best-known objects. 

The observations of nearby low-mass stars have been 
greatly improved, and their number has substantially in- 
creased, during the last few years. Distances have been 
determined with a very high precision with the Hippar- 
cos satellite. In parallel, ground-based measurements have 
provided high resolution spectra and multicolour photom- 
etry of an appreciable number of stars of the solar neigh- 
bourhood. On the other hand, models of stellar atmo- 
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spheres have largely benefited from progresses in the the- 
oretical description of microscopic physics, in particular 
opacities. The analysis of the observational data, using 
model atmospheres, has provided bolometric magnitudes, 
effective temperatures and abundances with an accuracy 
which had never been reached previously. In addition, sev- 
eral nearby stars are members of visual binary systems and 
their mass is known. 

The study of stars of low mass has many important 
physical and astrophysical implications. First, several un- 
certainties remain on the physics of their internal struc- 
ture. For instance, the unperfect understanding of the 
convective transport or of the non local thermodynami- 
cal equilibrium (non-LTE) effects in atmospheres and en- 
velopes have important consequences: atmospheres and 
envelopes serve as external boundary conditions for the 
interior models, fixing the radius, and are also used in the 
analysis of observational data to determine abundances, 
effective temperatures and gravities. On the other hand, 
transport processes at work in the deep interior need to 
be better constrained; they link the convective zone to 
deep nuclear burning regions, involving consequences for 
the surface abundances (lithium and light elements but 
also heavy elements) , and for the estimated age and global 
parameters, when fresh helium is brought to the burning 
core. 

An important question is also that of helium content, 
not measurable in the photosphere due to the lack of lines 
in the spectrum. The initial helium content of a star de- 
termines its lifetime and internal structure and is also a 
witness of past galactic history. The knowledge of the ini- 
tial helium abundance of stars born in different sites with 
different metallicities is therefore fundamental for studies 
of the chemical evolution of the Galaxy. 

We study here a sample of a hundred nearby disk stars, 
of spectral types in the range F to late K, by means of 
theoretical stellar models. Our aim is to discuss the abil- 
ity of the models to reproduce high quality observations 
and, when possible, to determine the helium content of 
the stars. 

In Sect. ^ we make a rapid summary of some results 
previously obtained on the subject, from studies of a few 
particularly well-known objects. In Sect. || we carefully 
examine the observational data available for nearby stars 
in order to extract subsamples corresponding to the very 
best accuracies. We then present the derived Hertzsprung- 
Russell (HR) diagrams of stars of the solar neighbour- 
hood which have the smallest error bars ever obtained. 
Section 4 describes the theoretical models used to com- 
pute isochrones for given chemical compositions. In Sect. 5 
we show evidence for discrepancies between observations 
and theoretical models, which are particularly obvious for 
stars with metallicities below [Fe/H] =-0.5. In Sect. 6 we 
examine this problem in the light of recent wo rks o n de- 
parture from LTE for iron (Thevenin & Idiart 199E ) and 
on microscopic diffusion of helium and heavier elements in 



metal-poor stars (Morel & Baglin 1999). We show that the 
cumulative effects of these two processes is large enough to 
remove the discrepancy shown in Sect. 5. Section 7 consid- 
ers the case of binary stars with known dynamical masses. 
Section 8 adresses the question of the helium content of 
individual objects of the sample. Section 9 summarizes our 
conclusions. 

2. Present status of helium content and 
mixing-length parameter determinations 

Since observations do not provide direct determinations 
of helium, assumptions have to be made for the initial 
helium content of models of low-mass stars. Very often 
it is supposed that the metallicity Z and helium Y are 
related through a constant AF/AZ -value. 

Ay/ AZ = (F - Yp )/Z is often determined from the 
results of the solar calibration and for any other star of 
known metallicity, the helium abundance is scaled on the 
solar one. In addition of the behaviour of Y with Z, an- 
other problem is the universality of the a = I /Hp ratio 
of the mixing length over the pressure scale-height in the 
mixing-length treatment of the convection. We now briefly 
review the various determinations of these parameters, in 
the closest stars and various astrophysical sites. 

2.1. Nearby stars 



Perrin et al. (1977) were the first to examine the HR di- 
agram of a selected sample of the nearest low-mass stars. 
Because in 1977 the error bars in the HR diagram were 
quite large, Perrin et al. (1977) suggested that all the non- 
evolved stars were sitting on the same ZAMS, whatever 
their metallicity was. They showed that this behaviour 
was reproduced by theoretical stellar models if helium and 
metallicity were related through a value of AY/AZ con- 
stant and equal to 5. Recently Fernandes et al. (1996) 



claimed that stars do not all lie on the same ZAMS and 
measured the observational lower main sequence width in 
the solar neighbourhood. They showed that if the width 
is entirely due to chemical composition dispersion in the 
solar neighbourhood, then AY/AZ is greater than 2 in 
the corresponding stars. 

More recently H0g et al. ( |1997D and Pagel & Portinari 
(|1998D have found AY/AZ = 3 ± 2 from a sample of 



nearby stars with Hipparcos parallaxes. 

We recall further constraints on these points. 

2.1.1. The Sun 

The Sun is a milestone in internal structure theory, be- 
cause its age and fundamental parameters are known with 
great accuracy, and helioseismology (Perez Hernandez & 



Christensen-Dalsgaard 1994) has brought additional con- 



straints. All observations have to be reproduced by the 
theoretical solar model, with underlying input physics and 
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free parameters. The initial helium abundance is a free 
parameter of the solar model and is estimated with the 



luminosity constraint (Christensen-Dalsgaard 1982). The 
difference between the present helium value derived from 
seismology(Y« 0.25) and the initial value obtained from 
calibration (Y w 0.275) provides a constraint on the input 
physics of the model: it can be explained by invoking the 
microscopic diffusion of helium and heavy elements which 
has taken place during the evol ution of the Sun (s ee e.g . 
Cox et al. |1989| , Richard et al. |l996| , Brun et al. 



1991) . 



Furthermore the mixing-length parameter, which enters 
the mixing-length phenomenology of convection is fixed 
by the radius constraint at a value « 1.7-1.8. 

2.1.2. Binary stars and clusters 

The study of a few visual binary systems of known mass, 
effective temperature, luminosity and metallicity provides 
further information on low mass stars. The two compo- 
nents of a binary system, which are assumed to have the 
same age and initial chemical composition, can be mod- 
elled simultaneously providing the age of the system, its 
initial helium content and the mixing-length ratio (Noels 
et al. |1991| ). The method was first applied to the a Cen- 
tauri s ystem by Noels et al. (1991). Recently, Fernandes 
et al. ( |l998[ ) studied four Population I low mass binary 
systems with high quality data and determined the initial 
helium abundance by mass for three of them with a preci- 
sion of 0.02. A global conclusion of these papers is that the 
mixing-length parameter a seems independent of metal- 
licity, and that a AY/ AZ ratio of 2.3 ±1.5 is appropriate 
(their Fig. 5). 

Using the high quality observations of Hipparcos to- 
gether with ground-based spectroscopic data. Ferryman 



et al. (1998) and Lebreton et al. (1997) determined the 
helium abundance by mass, Y, of the Hyades with a pre- 
cision of 0.02 . This abundance does not follow the over- 
abundance of the metallicity ( -I-.15 dex). This suggests 
that a scatter does exist on the AY/ AZ ratio. This paper 
also confirms the universality of the a parameter. 

2.2. Extragalactic observations 

The studies mentioned above have provided the initial he- 
lium content of a few stars of known metallicity which 
allows to calculate their helium to metal enrichment ra- 
tio AY/AZ = (Y -Yp )/Z. The primordial hehum abun- 
dance Yp has been determined by several groups: for in- 
stance, Balbes et al. ( |1993D g ave Y p =0.227±0.006 while, 
more recently, Izotov et al. (1997) found a higher value 
Yp =0.243±0.003. However, Hogan et al. ( [1997D consider 
this value as an upper limit. This gives two secure points 
of the Y(Z) relationship: Yp « 0.24 at Z = 0., and 
Yq = 0.275 at Z = Zq. The observations of blue compact 
galaxies and other systems have allowed to determine a 
AY/AZ ratio : a recent discussion is given by Izotov et 



al. (1997). The value is 2.3 ± 1.0, considerably below for- 
mer determinations « 4.0. 



2.3. Nucleosynthesis 

On the other hand nucleosynthetic predictions integrated 
over the whole stellar mass range lead to AY/AZ -values 
ranging from about 1 to about 2 depending on the in- 
clusion of stellar winds; AY/AZ can even reach values 
of about 5 if black holes are considered (Maeder 1992). 
Although helium is expected to increase with metallicity, 
it must be noted that metals are only produced by SNe, 
whereas helium is also produced by mass-loss of medium- 
mass stars, in their post-AGB phase. 

This body of results indicates that a significant dis- 
persion of the helium abundance around the solar-scaled 
value cannot be excluded. 



3. Observational determination of the stellar 
fundamental parameters 

We study an homogeneous sample of 1 14 late-type nearby 
stars of spectral types in the range from F to late K, care- 
fully selected by M.-N. Perrin. They are part of the pro- 
posal 132 (M.-N. Perrin), accepted as an Hipparcos pro- 
gram in 1982. Later on, the proposal was updated as the 
INCA proposal Oil by A. Baghn, M.-N. Perrin, Y. Le- 
breton and R. and G. Cayrel. These stars are closer than 
about 25 parsecs which ensures an excellent accuracy of 
their parallax determination by Hipparcos. 

Among these, we have retained the stars which have 
been submitted to detailed spectroscopic analysis from the 
ground and which appear in the last version of the Cata- 
logue of [Fe/H] determinations: 1996 editionhy G. Cayrel 
de Strobel et al. ( 1997| ). Their metal to hydrogen ratio 
[Fe/H], i.e. the logarithm of the iron to hydrogen ratio 
(by number of atoms) relative to the solar value, ranges 
from about -1.0 to about -1-0.3 dex corresponding to Pop- 
ulation I and thick-disk population. For each star we as- 
sign an averaged abundance determined from the above 
catalogue (for more details see the caption of table 1). 
The mean internal error on [Fe/H] is of the order of 0.1 
dex . This does not include systematic errors, as the fact 
that all abundances in the catalogue are not corrected for 
NLTE effects. In determinations based on spectra taken 
with solid state detectors, with high quantum efficiency, 
the random error coming from equivalent width measure- 
ments enters for about 0.03 dex only, the remaining part 
coming from uncertainties in the fundamental parameters 
Tcff and \ogg. 

We eliminated the suspected unresolved binaries and 
only kept stars with parallaxes determined with an accu- 
racy better than 5 per cent. Then, among the remaining 
stars, we extracted three homogeneous and independent 
subsamples. 
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Table 1. Observational parameters for Samplel. Parallax tt, and relative error on it CTtt/t'' come from the Hipparcos 
main catalogue. The apparent magnitude V is from the Hipparcos Input Catalogue. Mboi and Toff were derived from 



the bolometric fluxes of Alonso et al. (1995, 1996a). Note that the zero-point of bolometric magnitudes is not that used 
by Alonso, but is such that Mi,oi =4.75 for the Sun. [Fe/H] values are a weighted average of individual values from 
spectroscopic analyses taken in Cayrel de Strobel et al. (1997). A quality index "qlt" was attributed to the adopted 
value of [Fe/H] according to the following code: 

4: average of at least six determinations obtained with recent high S/N spectra 
3: average of at least three determinations obtained with recent high S/N spectra 
2: at least one determination obtained with recent high S/N spectra 

1: [Fe/H] based on photographic spectra, keeping only high quality work, shown to be exempt of large systematic 
errors (often from Hernshaw, see Fuhrmann 199 j ) 
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Sample 1 is constituted of 33 stars with directly de- 
termined bolometric fluxes, so the problem of the value 
of the bolometric correction is eliminated. The relevant 
data for Sample 1 are li sted i n Table [l]. The fluxes were 
derived by Alonso et al. (1995) with an accuracy of about 
2 per cent by integrating UBVRIJHK photometry. The 
average absolute error on bolometric magnitude resulting 
from the parallax and bolometric flux uncertainty is of 
about 0.03 to 0.06 magnitude. The effective temperature 



was obtained by Alonso et al. ( 1996a ) from the bolometric 
flux by the I nfra-Red Flux Method (IRFM) (Blackwell & 
Shallis 1977 ) . The method makes use of a grid of theoreti- 
cal model line-blanketed flux distributions (Kurucz 1991). 
The mean internal error on effective temperature is about 
80 K f or temperatures greater than 4200 K (Alonso et al. 
1996a ). Fig. l| shows the resulting positions of the stars 
in the (log T^g , Mboi ) diagram with their individual er- 
ror bars. The striking feature of Fig. |l| is the lack of clear 
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Fig. 1. The Hipparcos HR diagram of the 33 best known 
nearby stars (Sample 1, see text). The accuracy on the 
parallax is better than 5 per cent and bolometric fluxes are 
available for each star. Individual errors bars are given for 
each object by the horizontal and vertical sizes of the cross. 
Filled circles, open circles and open triangles respectively 
correspond to stars with [Fe/H] values in the interval [- 
0.10, 0.30], [-0.45, -0.10] and [-1.05,-0.45] 



separation between stars having a solar metallicity, and 
those in the metallicity range [-1.0;-0.5], confirming with 
high quality data the result already obtained by Perrin 
et al. ( 1977| ), on less accurate data. In Fig. |^ the sample 
of stars plotted in Fig. |l] is split into 2 subsamples corre- 
sponding 1) to stars of solar metallicity, plotted in Fig. 2a 
and 2) to moderately deficient stars, plotted in Fig. 2b. 
Theoretical isochrones are also plotted in Fig. |^ but this 
will be discussed later in Sect. ^. 

Sample 2 is constituted of 64 stars with effective tem- 
peratures derived from detailed spectroscopic analysis. For 
each star at least one determination of the effective tem- 
perature and of [Fe/H] was made with modern detectors 
(CCD or Reticon). We adopt a mean internal uncertainty 
on effective temperature of 150 K. Bolometric magnitudes 
were obtained using bolometric corrections from Alonso 
et al. (1996a, 1996b| ), with a correction of the zero point 
giving a bolometric magnitude of 4.75 to the Sun. The 
resulting positions of the stars in the (log T^s , Mboi ) 



diagram are plotted in Fig. ^ and will be discussed later 
in Sect. |. 

Sample 3 is constituted of 15 stars with effective tem- 
peratures derived from the {V — K) color index an d from 
[Fe/H] accordi ng to the results f rom C arney et al. (1994), 
Johnson et al. (|l966|) , Koorneef (|l983|) . The mean internal 
uncertainty on effective temperature is of about 75 K. The 
bolometric magnitudes are derived from the V magnitude 
given in the Hipparcos input catalogue, using bolomet- 
ric corrections from Alonso et al. ( |l9964 |l996b|) , with the 
same zero-point. Fig. ^ shows the resulting positions of the 
stars in the (log T^g , Mboi ) diagram with their error bars 
and will be discussed in Sect. |. Very clearly Fig. ^ has 
a much larger scatter than Fig. l| 2 and 4, showing that 
the effective temperatures derived from detailed analyses 
are not tightly connected to the true effective tempera- 
tures. The effective temperatures derived from {V — K) 
on the contrary are tightly bound to those obtained by 
the IRFM, as it can be directly checked in the full Alonso 
sample. It is not surprising, as effective temperatures from 
detailed analyses are affected by other parameters: grav- 
ity, departures from LTE in ionisation equilibria, and a 
variety of different techniques in fixing the triad (Tcff , log 
g and [Fe/H]). 

4. Theoretical stellar models 

J^.l. Input physics 

The stellar evolution calculations were computed with the 
CES AM code (Morel |1997| ) in which we have included ap- 
propriate input physics for the mass range under investi- 
gation. 

The CEFF equatio n of state (Eggleton et al 
Christensen-Dalsgaard 1991), was used. It includes the 



1973 



Coulomb corrections to the pressure and is appropriate 
when modeling low-mass stars of m ass greater than about 
0.6 Mq (Lebreton & Dappen 



We used the nuclear reaction rates given by Caughlan 
& Fowler ( p^88| ). 

We calculated the initial compo sition of the models 
either from the Grevesse & Noels ( 1993 ) solar mixture 
(GN93 mixture) or from a GN93 mixture where the a- 
elements O, Mg, Si, S, Ca, Ti are enriched relative to 
the Sun ([a/Fe]= -1-0.4 dex). Such an enrichment of a- 
elements is observed in metal deficient stars with metallici- 



ties [F e/H] lower than -0.5 (Wheeler et al. |l989| , Fuhrmann 
1998) and its effect on models has to be taken into account. 

We took the most recent OPAL opacities (Iglesias 
& Rogers 1996) provided for the two mixtures consid- 



ered. They were complemented at low temperatures by 
atomic and molecular opacities from Alexander & Fergu- 
son ( |1994D for the GN93 mixture or from Kurucz ( |199lD 
for the a-enriched mixture. Low and high temperatures 
tables were fitted carefully at a temperature of about 
10 000 K, depending upon the chemical composition. 
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In most models the atmosphere was described with 
an Eddington's T{t) law which is easy to use. In order 
to estimate the uncertainties resulting from this choice 
we built a few models using T(t) laws derived from the 
ATLAS9 atmosphere models (Kurucz |199l| ) . This requires 
to calculate detailed atmosphere models with ATLAS9 for 
many values of the gravity, of the effective temperature 
and of the metallicity and then to derive the corresponding 
T{t) laws (see Morel et al. 1994). Interpolation of the 



Y{z) = yp + (z/Zq) {Yq - yp) 



19931) . This implies a 



T(t) laws is then performed to model the considered star. 
We checked on a test-case, that the use of those better 
atmospheres does not shift the position of the model in the 
HR diagram by an amount larger than the observational 
error bars. 

Convection was treated according to the mixing-length 



theory of Bohm-Vitense (1958). The mixing- length pa- 
rameter Of, ratio of mixing- length to pressure scale-height, 
is a free parameter of the models. As shown by Fernan- 
des et al. (1998) the Sun and four visual binary systems 



spanning a wide range of masses and metallicities could be 
calibrated with very close values of a. Moreover the slope 
of the Hyades main-sequence i s well- reproduced with a so- 
lar a- value (Fe rryma n et al., 1998). On the other hand, 
Ludwig et al. ( |l997 ) made a calibration of the mixing- 
length based on 2-D radiation hydrodynamics simulations 
of solar-type surface convection and found that variations 
of a of about 10 per cent around the solar value are ex- 
pected in the domain of effective temperature, gravity and 
metallicity studied here. We adopt the solar mixing-length 
value in our calculations. For unevolved stars a change of 
£/Hp of ± 0.15 produces a change in effective temperature 
of the order of 50 K, smaller than the mean observational 
error on Tog. 

With the input physics described above the calibra- 
tion in luminosity and radius of a solar model having 
{Z/X)q =0.0245, where X is th e hydrogen abundance by 
mass (Grevesse and Noels 1993) requires a mixing- length 



£ = 1.64 Hp, an initial helium abundance y=0.266 and a 
metallicity Z0=O.O175. 

4-2. Grids of stellar models and isochrones 

In order to determine the ZAMS position , we calculated 
zero and terminal age main sequences (ZAMS and TAMS) 
for 10 mass values ranging from 0.5 to 1.4 Mq , 5 helium 
values ranging from 0.18 to 0.43 and metallicities values 
Z= 0.004, 0.007, 0.01, 0.015, 0.02, 0.03, 0.04 and 0.06. 

To discuss the global features of the HR diagram and 
the position of particular objects, we calculated detailed 
evolutionary sequences from the ZAMS to the beginning 
of the red-giant branch, for 14 masses ranging from 0.5 to 5 
Mq and we derived isochrones using the Geneva isochrone 
program (Maeder 1974). We chose metallicities corre- 
sponding to the observational range (i.e. [Fe/H]=-|-0.3, 
0.0, -0.5 and -1.0) and solar-scaled values of the helium 
abundance given by : 



with yp = 0.227 (Balbes et al. 
AY/AZ value of 2.2. 

Two distinct grids of models were calculated in the 
metal deficient cases ([Fe/H] = -0.5 and -1.0): a grid with 
normal solar mixture and a grid with an a-elements en- 
hancement of 0.4 dex. For the solar mixture grid the metal- 
licity Z in mass fraction is related to [Fe/H] by: 

[Fc/H] =log(Z/X)-log(Z/X)0 

where {Z/X)q is the ratio of the solar mixture considered. 
For the a -elements enhanced mixture the relation be- 



[Fe/H] = log{Z/X) - log iZ/X)Q - 0.305 

where 0.305 = log {Z/X)q^^ - log {Z/X)q, difference be- 
tween a solar a-elements enhanced mixture and a 'normal' 
solar mixture. 

The detailed results of all the stellar models computed 
for that work will be published separately (Lebreton et al. , 
in preparation) and are available on request. 

5. Comparison between theoretical isochrones 
and observational data 

We present in Fig. 2a and in Fig. 2b the HR diagram of 
the 33 stars of Sample 1. We have split the sample in two 
subsamples: Fig. 2a corresponds to stars of solar metal- 
licity or close to it (-0.45 < [Fe/H] < + 0.30) while 
Fig. 2b is constituted of the moderately metal deficient 
stars (—1.05 < [Fe/H] < —0.45). We have superimposed 
on those diagrams theoretical isochrones associated to the 
limits of the metallicity ranges considered. Isochrones with 
[Fe/H] values of -0.5 and -1.0 have been derived from mod- 
els calculated with an a-elements enriched mixture and 
have solar-scaled helium values, y=0.256 and 0.236 re- 
spectively. The metal-rich isochrone has a solar a-elements 
ratio and y=0.32. We have also plotted the position of the 
solar ZAMS on Fig. 2a. Isochrones are given for ages in 
the range 8-10 Gyr representative of the age of the galactic 
disk. 

Fig 2a shows that the stars of solar metallicity or close 
to it lie in the region defined by the theoretical models 
corresponding to their associated metallicity range. We 
notice the particular position of 7 Lep (HIC 27072) which 
is a very young star lying close to the ZAMS. 

However, it is worth to note that the slope of the the- 
oretical main sequence agrees well with the observational 
slope. This, once more, is in favor of the uniqueness of the 
mixing-length parameter in low-mass stars. Furthermore 
the width of the global observed main-sequence band for 
metal-rich stars, which is around 0.3 mag, and the theo- 
retical width corresponding to the same metallicity range 
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Fig. 2. Same as Fig. |l| (i.e. stars of Sample 1) but split in two metallicity domains: stars of solar metallicity and close 
to it (Fig. 2a) and moderately deficient stars (Fig. 2b) . Theoretical isochrones are overlaid on the observational data. 
Fig. 2a: the lower isochrone (8 Gyr) is for [Fe/H]=— 0.5, y=0.256 and [a/Fe]=0.4; the upper isochrone (10 Gyr) is for 
[Fe/H]=0.3, Y—0.32 and [a/Fe]=0.0. The dashed line is a solar ZAMS. The young star 7 Lep is the brightest star in 
the figure. Fig. 2b: the lower isochrone (10 Gyr) is for [Fc/H]=— 1.0, y=0.236 and [a/Fe]=0.4; the upper isochrone 
(10 Gyr) is for [Fc/H]=— 0.5, F=0.256, [a/Fc]=0.4. All stars, but one, are lying above the lane defined by these two 
isochrones. 



with solar-scaled helium abundances are quite similar (see 
also Fig ^ . 

Turning now to the moderately deficient stars plotted 
in Fig. 2b we note that all the stars but one are outside the 
theoretical band corresponding to their metallicity. Stars 
are located on theoretical isochrones with a metallicity 
higher than their observed metallicity. 

The same behaviour is seen in the two other selected 
samples in which the effective temperatures of stars have 
been determined independently by different authors (see 
Sect. H). Fig. 3a and 3b are similar to Fig. 2a and 2b but for 
the 64 stars with effective temperatures derived from de- 
tailed spectroscopic analysis (Sample 2) while Fig. 4a and 
4b represent the situation of the 15 stars of Sample 3. Al- 
though the error bars on effective temperature are greater 
in Fig. 3a and 3b and although they are few objects in 
Fig. 4a and 4b, the same tendency is found, strengthen- 
ing the conclusion that classical theoretical isochrones and 



actual observations of moderately metal deficient stars do 
not match. 

Note that, in the confrontation of observational posi- 
tion of stars with stellar isochrones, errors on [Fe/H] and 
[a/Fe] which are in the range 0.10-0.15 dex, are respon- 
sible for a hidden enlargment of the error boxes, because 
the chemical composition of the models to be compared 
with the observed stars are not known exactly. 



6. Resolving the metal-poor discrepancy 

The most natural attempt to solve the discrepancy is by 
trying to adjust the "free" parameters, the helium con- 
tent Y being the first target. However, for the stars in 
the metallicity bin [—1.05, —0.45], a Y < 0.2 would be re- 
quired, a cosmologically unacceptable value. A change in 
a = (/Hp does not help either, as it is producing a neg- 
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Fig. 3. Same as Fig. || for stars of Sample 2 (Mboi is obtained from the Hipparcos parallax and V magnitude and Toff 
from detailed spectroscopic analyses, see Sect. O) 



ligible effect in the low part of the main sequence under 
consideration here. 



So, as noted in Lebreton et al. (1997) the discrepancy 
in the expected and actual positions of metal-poor stars 
needs either a drastic change in the zero-point of effective 
temperature (of the order of 200 K) or accepting the view 
that other processes must be included in the interpretation 
of the observations. It would be very difficult to plead a 
large correction in the effective temperatures at metallicity 
[Fe/H] = —0.7 and none at solar metallicity. So, we do not 
further consider this explanation and, now focus on the 
second one. 

Two processes are considered in this section: (i) de- 
partures to LTE in the determination of [Fe/H] (ii) micro- 
scopic diffusion of helium and other heavy elements. 

6.1. non-LTE departures for iron 

All the [Fe/H] values we have used were computed with 
the LTE approximation. Computing NLTE abundances 
for iron is a formidable task. However, it must be noted 



that Bikmaev et al. ( 1990 ) have claimed that a substantial 
overionisation occurs for Fe I in subdwarfs, and only Fe H 



lines should be used for determining iron abundances. 
Holweger and coworkers (private communication, 1993) 
found the same phenomenon in trying to compute depar- 
tures from LTE for the very m etal-p oor star CD — 38°245. 
More recently, Nissen et al. (1997) found that spectro- 
scopic gravities in subdwarfs were systematically lower 
that the gravities derived from Hipparcos distances. Fi- 
nally, Thevenin & Idiart ( 199!: ) have undertaken NLTE 
computations of iron abundances in one hundred metal- 
poor stars, and found non-LTE corrections for several of 
our stars, of the order of 0.15 dex, for the mean metallic- 
ity ([Fe/H] = -0.72) of our -1.05 < [Fe/H] < -0.45 
subsample. Non-LTE corrections are negligible for stars 
with solar metallicities (same reference). The metallicity 
of isochrones to be compared to the set of stars with a 
mean LTE metallicity [Fe/H] =-0.72 is therefore the true 
NLTE metalhcity [Fe/H]=-0.57. The situation is improved 
but a clear departure is still existing (Fig. ||). 

6.2. Sedimentation of heavy elements 

The effects of sedimentation of heavy ele ment s have re- 
cently been computed by Morel & Baghn (1999), for stars 
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Fig. 4. Same as Fig. ||for stars of Sample 3 (Toff is obtained from the (V-K) color index and from [Fe/H], see Sect. ^ 



of metallicity [Fe/H] = -0.7, -0.9, -1.2,-1.7, and for 

5 stellar masses, 0.85, 0.80, 0.75, 0.70, and 0.60 Mq. To 
produce significant effects in the HR diagram microscopic 
diffusion has to proceed on a sufficiently long timescale. 
This condition is fulfilled for the thick disk subsample un- 
der consideration here: the 10 Gyr isochrones of the Morel 

6 Baglin paper are a suitable choice for our subsample. 
Two effects are present: 

After a time of 10 Gyr, the stratification of the heavy 
elements translates the evolution point in both logTeff 
and Mboi- This effect is small: at 0.7 Mq and an ini- 
tial metallicity of [Fe/H] ~ —0.75 the point is moved by 
AMboi = -0.029 and A log Toft = -0.0004. But the sur- 
face metalhcity of the star has dropped to [Fe/H] ~ —0.81. 
Therefore, the theoretical position of a star of present 
metallicity [Fe/H] = —0.75, must be computed with a 
modified initial composition, leading after a 10 Gyr evolu- 
tion, to the metallicity we observe now. This is the largest 
part of the correction, called the calibration correction 
by the authors. From their table 2, and their Fig. 5 we 
have computed the corrections to be applied to a stan- 
dard isochrone for the mean metallicity ([Fe/H] = —0.72) 
of our sample in the range —1.05 < [Fe/H] < — 0.45. 



6.3. Cumulated effect 

Applying these corrections to the dashed curve in Fig. |^ 
produces the dot-dashed isochrone. No systematic depar- 
ture appears any more between the observations and the 
theoretical isochrone corresponding to the present surface 
metallicity of the stars. The initial discrepancy noted at 
the beginning of this section is then completely removed. 
The only star presenting a 2a discrepancy is the star 
HD 132142 at (3.707, 5.697). But there is a single spectro- 
scopic determination of its metallicity and its photometry 
indicates a solar metallicity. So this exception is not a se- 
rious worrying. 



7. Binary stars 

One more constraint can be obtained from the stars of the 
sample which are component of a binary system having 
known orbit and masses. One of us (JF) has already used 
a few of the best known binaries in order to check the pre- 
dictions of theoretical evolution computations (Fernandes 



et al. 1998). It is worth to see if, when known, the dynam- 



ical mass of the object confirms, or disproves, its location 
along the proper isochrone obtained above. The most in- 



10 



Y. Lebreton et al.: The HR diagram for late- type nearby stars after Hipparcos 




3.75 



3.65 



bgTeff 



Fig. 5. This figure is a large scale HR diagram for the sub- 
sample of unevolved stars in the LTE metallicity range - 
1.05 < [Fe/H] < -0.45. The full hne is a standard isochrone 
for the mean LTE metallicity of the sample [Fe/H] = -0.72. 
The dashed line is the standard isochrone corresponding to 
the NLTE mean metallicity of the sample [Fe/H] — -0.57. 
The dot-dashed line corresponds to an isochrone includ- 
ing microscopic diffusion of the elements for an age of 10 
Gyr. The surface metallicity is the NLTE value [Fe/H]=- 
0.57, but the initial metallicity, which is very close to the 
mean interior metallicity was [Fe/H] w —0.5. The fit is 
now satisfactory. The tick marks show the value of the 
mass along the upper isochrone. The tick mark near 85 
Peg A corresponds to 0.85 Mq (not labeled for clarity) 



teresting case is the binary /i Gas, the only metal-poor 
star, having a well determined mass. 



7.1. /i Cas 



Drummond et al. (1995) have derived the mass of /i Gas 
from its orbit, with the help of speckle interferometry to 
see its very faint companion (Am w 5.0). The mass of 
fi Gas A, readjusted with the better distance provided 
by Hip parcos is 0.757 ± 0.06 Mq. The Alonso et al. 
( |1996aD effective temperature for (A-fB) is 5315 ± 82 K, 
which corrected for the presence of the companion be- 
comes 5339 ± 85 K, in good agreement with Fuhrmann 
( |1998D effective temperature of 5387 ± 80 K. We obtain 



Mboi = 5.634 ± 0.035 for /i Gas A from its parallax, and 
again correcting for the presence of the companion. The 
position of the star is plotted with respect to the uncor- 
rected and corrected isochrones in Fig. ^. The evolution- 
ary masses along the dot-dashed isochrone are shown with 
tick-marks. The location of fi Gas is near 0.8 Mq, within 
the limits of the dynamical mass 0.757 ± 0.06 Mq. Un- 
fortunately, the error bar on the mass is still a bit large 
to add an interesting new constraint. The influence of the 
uncertainties on the values of the various parameters of 
the models on the position of /i Gas in the HR diagram 
are documented in table 2 (see Lebreton 199S for more 
details). 



Table 2. budget of the impact of parameter variations on 
the position of a star like n Gas in the HR diagram. The 
reference model has an age of 12 Gyr, a mass of 0.757 Mq, 
Y = 0.245, [Fe/H] = -0.86, and [a/Fe]= 0.3 dex 



parameter 


increment 




SMboi 


mass(MQ ) 


0.06 


0.031 


-0.607 


age (Gyr) 


2 


0.005 


-0.150 


Y 


-0.015 


-0.012 


0.172 


[Fe/H] (dex) 


-0.15 


-0.015 


0.188 


[a/Fe] (dex) 


0.10 


-0.004 


0.064 


a (MLT) 


-0.15 


-0.004 


0.004 


atmosphere 


Eddington/Kurucz 


-0.001 


0.001 


diffusion 


yes/no 


-0.007 


0.040 


low-T opacity 


-f 30 per cent 


-0.049 


0.742 



7.2. 85 Peg 
Fernandes et al 



showed that, using the data avail- 
able on this object, there was no theoretical solution fittin g 
the two components . Later on, Martin & Mignard (1998), 
have restudied several binaries from Hipparcos data, and 
shown that it is very difficult to escape the conclusion that 
the mass of the primary and the secondary are very sim- 
ilar, notwithstanding the fact that the secondary is 3.2 
magnitude fainter than the primary. New determinations 
of the effective temperature of 85 Peg A by Thevenin & 
Idiart (1999) and G. Van't Veer (private communication, 
1998) confirms that 85 Peg A has a metallicity very simi- 
lar to that of Gas A, an effective temperature of 5550 ± 
100 K, and Mboi=5.22 ± 0.06. Its location in the HR dia- 
gram (Fig. ^) is pratically on the same isochrone as ^ Gas, 
and corresponds to a mass of 0.85 Mq, higher than the 



mass given by Martin & Mignard ( |1998D (0.705 Mq), but 
below the mass 0.95 ± 0.1 given by Duquennoy & Mayor 
( 1991[ ) from the spectroscopic orbit. 85 Peg A does not 
deviate any more from the proper isochrone, once its high 
[a/Fe] ratio, its true iron abundance corrected for NLTE 
effects, and sedimentation of heavy elements, are all taken 
into account. The problem is with 85 Peg B, which is too 
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massive for its absolute magnitude. The colour {V — K) 
of 85 Peg B can be derived from the V and K magnitudes 
of (A+B) , V=7.75, K==3.94 (Johnson et al. [1968D , the 
magnitudes Va= 5.81 and Vb= 9.0, and assuming that 
iy — K)a = 1.59 ±0.02 from its effective temperature and 
Alonso calibration of {V — K) versus (Tf-fj, [Fe/H]). The 
result is {V — K)b— 3.46 ±0.06, corresponding to about 
Tf.ff = 3950 K. The associated mass and absolute bolo- 
metric magnitude are respectively 0.54 (smaller than the 
dynamical mass of about 0.7 to 0.8 Mq ) and 8.0. With a 
bolometric correction of -1.1 and the transformation from 
apparent to absolute magnitudes this gives V_b — 9.57 too 
faint by 0.6 magnitude. One possibility is that 85 Peg B 
is a spectroscopic binary itself, but we shall leave this hot 
subject for another paper, this one having at least clarified 
the case of 85 Peg A. 

7.3. T] Cas 

Nothing really new, not already in Fernandes et al. ( 1998D , 
can be said for this star. The fact that it is more evolved in- 
troduces one more parameter (age) and prevents to check 
the unevolved position of the star. The star is neverthe- 
less identified in Fig. and does not raise a particular 
problem. 

8. How does helium content vary v^rith 
metallicity? 

Our isochrones are computed with the assumption that 
the helium content varies in proportion with metallicity, 
namely: 



Y{Z) = Fp + {Z/Zq) {Yq - Fp) 



(1) 



Any unevolved star presenting a significant deviation 
with respect to the isochrone corresponding to its metal- 
licity may indicate that its helium content deviates from 
this simple, naive assumption. Actually there is already an 
indication that t he Hyades do not follow this rule (Perry- 
man et al. 1998| ), Hyades stars being metal-rich by 0.15 
dex with respect to the Sun, without being correspond- 
ingly helium-rich. On the contrary, the metal-rich stars of 
the a Cen system were found also helium-rich (Y=0.30) by 
Fernandes & Neuforge, (1995) in agreement with equation 



(1). 

In Fig. H we have plotted the most unevolved stars 
(Mboi > 5.2), labeled by their value of [Fe/H], with a solar 
ZAMS and a grid of isochrones of 8 Gyr, of metallicities 
—0.3, 0. and +0.3. Only stars of metallicity corresponding 
to this interval are plotted. We admit that the age of these 
thin disk stars is between and 8 Gyr. So if the error 
bars intersect either the right metallicity ZAMS or a right 
metallicity isochrone there is no compelling evidence that 
the scaled helium value is wrong. Globally all stars within 
the metallicity range -0.3 < [Fe/H] < 0.3 are within the 
corresponding isochrones. But inside the lane they define. 
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Fig. 6. This figure is a large scale HR diagram for the sub- 
sample of unevolved stars in the LTE metallicity range - 
0.3 < [Fe/H] < 0.3. Each star is labeled by its LTE [Fe/H] 
value, followed by the HD number of the star. The full 
line is a standard isochrone for the metallicity [Fe/H] — 
0. The dashed line is the standard isochrone correspond- 
ing to the LTE metaUicity of the sample [Fe/H] = -0.3 
The dot-dashed line corresponds to an isochrone of metal- 
licity +0.3. All isochrones are for an age of 8 Gyr , but 
evolutionary effects are very small, and age is almost ir- 
relevant. If stars stand between the limit isochrones, they 
do not always closely match the isochrone corresponding 
to their metallicity, leaving the impression that AY/AZ 
is not constant in the sample. 



some stars are not located exactly where one would expect 
to see them. We note that 5 stars (HD 37394, HD 131977, 
HD 166620, HD 201091 A, HD 219134) with metaUicity 
very close to zero or even slightly negative, stand on the 
[Fe/H] = +0.3, 8 Gyr isochrone. The non-LTE correction 
for their iron abundance is negligible at this metallicity. 

As disk stars are younger than halo ones, sedimenta- 
tion is less important. Its effect has been estimated, using 
a 5 Gyr, 0.8 Mq star of solar composition as a template. In 
this star, the shift in the HR diagram due to microscopic 
diffusion is of approximately 50K in effective temperature 
and 0.02 dex in luminosity. These values are comparable 
to the observational error bars on these parameters for 
individual objects. Only for the oldest and more massive 
stars, the shift can reach more important values, around 
120 K. Unfortunately, ages are most of the time unknown 
for individual unevolved stars of the disk. 

Though certainly responsible for a scatter in the HR 
diagram (due to mass and age differences), and for a sys- 
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tematic very small shift towards lower effective tempera- 
ture, any clear signature of this effect in the data seems 
impossible to decifer presently. Further work is planned 
on this point. 

Young stars should not show this effect. Indeed, in an- 
other paper (Pcrrymann et al.l998 , the content in helium 
of the Hyades, sufficiently young to forget about diffusion, 
has been determined, anf found to be solar like or slightly 
below, but certainly not enhanced at the level of the .15 
dex level of [Fe/H]. Also the very young active stars as 
e Eridani (HD 22049) or HD 17925 do not significantly 
deviate from the solar ZAMS. 

Another possibility for the 5 stars close to the -f .3 
dex isochrone is that these stars have less helium than the 
scaled value given by equation (1). Interpolation of ZAMSs 
corresponding to various (Y,Z) combinations could pro- 
vide an estimate of Y, for a star of known Z. However 
we must remember that the error bar of 0.1 dex on [Fe/H] 
translates into an error of 0.015 in Y, because of the (Y,Z) 
degeneracy in the HR diagram. Therefore we have decided 
to refrain from giving an helium content for individual un- 
evolved stars until the effects of sedimentation are prop- 
erly quantified at high metallicities. 

Summarizing, our data are consistent with the scaled 
helium content relationship at a coarse level of accuracy, 
of the order of a factor of two in metallicity, but suggest 
that a scatter exists at a finer resolution. 



9. Conclusions 

The location in the (log Tea, Mboi) plane of a sample 
of over one hundred nearby stars, covering a metallicity 
range of —1.0 < [Fe/H] < -|- 0.3, with Hipparcos 
parallaxes of relative accuracy better than 5 per cent, has 
been obtained. Most stars have bolometric magnitudes di- 
rectly measured, and effective temperatures derived from 



the Infra-Red Flux Method, by Alonso et al. ( |1996aD . The 
metallicities (LTE) have been derived from the 1996 edi- 
tion of the catalogue of [Fe/H] determinations by Cayrel 
de Strobel et al. ( |1997| ). 

The HR diagram shows an unexpected clumping of 
stars of low metallicity, only slightly separated from the 
solar metallicity sequence, as found much earlier, from 



ground-based parallaxes by Perrin et al. (1977). However 
the explanation of this fact by a variation of the helium 
content with metallicity Ay/AZ = 5.5 is no more ten- 
able, with the narrow range of variation of Y between 
Yp Ki 0.24 and Yq 0.28. A comparison of the data 
with the standard theoretical isochrones, computed with 
improved opacities, show that unevolved stars with solar 
metallicities agree more or less with the theoretical ex- 
pectations, whereas stars with metallicities below -0.5 dex 
deviate from the isochrones having the same metallicity. 
The case of ^ Cas, a binary of metallicity —0.7, is partic- 
ularly illustrative. 



Two effects, recently studied in the literature have 
been investigated to explain this observation. The first 
one concerns only the atmosphere of the star. Until now, 
non-LTE abundances for iron were not currently available 
in the literature. Thevenin & Idiart ( 1999| ) have studied 
more than one hundred stars, including /i Cas and many 
other metal-poor stars, and found a non-LTE correction of 
+0.15 dex for [Fe/H] in /x Cas. Morel & Baglin (|1999D have 



studied the microscopic diffusion of helium and heavier el- 
ements in metal-poor stars of various mass. They found 
that the atmosphere is depleted in iron after 10 Gyr, and 
that the effective temperature of the star is shifted by 
about —30 to —200 K for its luminosity. The combination 
of the two effects brings Cas back on its theoretical se- 
quence. The bulk of the combined correction is due to the 
fact that both the non-LTE correction, and the diffusion 
of heavy elements, make the atmospheric LTE iron abun- 
dance not representative of the mean inside abundance of 
iron. 

Among stars with metallicities closer to the solar 
metallicity there is a global agreement between the ob- 
servations and the isochrones with a scaled helium abun- 
dance, although some scatter may be present. New com- 
putations of sedimentation at these metallicities would be 
necessary to further resolve this issue. 

Also an enlarged sample of stars, with distances known 
with an accuracy of 1 or 2 per cent, appears necessary for 
a good determination of the helium content of individual 
stars. The GAIA mission is very promising in this respect 



(Ferryman et al. 1997). With GAIA the mean precision 



on parallax measurements will be at the 10 micro-arc- 
second (mas) level up to V=15, to be compared with 1 
mas with Hipparcos up to V=9. In parallel, progresses in 
the analyses of stellar spectra are relevant, as exemplified 
by the introduction of departures from LTE of the iron 
atom. 

More stellar masses are also deeply needed, not a single 
mass being known for stars more metal-poor than fi Cas. 
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